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Résumé 

Dans la première partie de cet article, on établit une liaison étroite entre la 
formule de Euler - Maclaurin et l'équation fonctionelle de Rota - Baxter ; ces deux 
choses étant plus ou moins équivalentes. 

Dans la deuxième partie, on présente une simple démonstration d'une formule de 
Ramanujan sur la sommation de certaines séries exponentielles. Ceci a fait l'objet 
d'un exposé à Toulouse, en mai 2007. 
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Première Partie 

ÉQUATION DE ROTA - BAXTER 
ET FORMULE SOMMATOIRE D'EULER - MACLAURIN 



§1. Définition de Jacob BernouUi 



1.1. Les nombres qu'A, de Moivre, puis Euler, ont appelés nombres de BernouUi, 
ont étés introduits par Jacob I BernouUi (1655 - 1705), dans son livre Ars Con- 
jectandi sur les probabilités, cf. [B], Pars secunda, Caput III, pp. 97 - 98. Ce livre 
a été publié en 1713, quand Euler avait 6 ans (Euler fut un élève du frère de Jacob, 
Johann, et un ami de ses deux fils, Nicolas et Daniel). 

BernouUi commence par un calcul de polynômes qu'il désigne par J rf ; nous 
convenons de la notation 

S^in) = V + + ... + n'' 

La méthode de calcul est basée sur le triangle de Pascal (qui à l'époque a servi 
pour la définition des numerorum figuratorum, alias coefficients binomiaux). Cette 
méthode était déjà connue de Pierre de Fermât. 

Voici ce qu'écrit BernouUi : 



"... Atque si porrô ad altiores gradatim potestates pergere, levique negotio 
sequentem adornare laterculum licet : 

Summae Potestatum 
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/n^o = ^nii + inio + f - 1 + 1 - in^ + A 



n 



66' 

Quin imô qui legem progressionis inibi attentuis ensperexit, eundem etiam con- 
tinuare poterit absque his ratiociniorum ambabimus : Sumtâ enim c pro potestatis 
cujuslibet exponente, fit summa omnium seu 



/ 



c • c — 1 • c — 2 • c — 3 • c — 4 ^ „ r 
2-3-4-5-6 

c-c— 1-c — 2 c — 3-c — 4-c— 5-c — 6^ ._7 . . 
+ 2.3.4.5.6.7.8 " • • • & demceps, 

exponentem potestatis ipsius n continue minuendo binario, quosque perveniatur ad 
n vel nn. Literae capitales A, B,C,D k. c. ordine dénotant coëfficientes ultimorum 
terminorum pro / nn, / n'^,Jn^,Jn^, & c. nempe 

A=-, B = - — , C=—, D = - — . 
6 30 42' 30 

Sunt autcm hi coëfficientes ita comparati, ut singuli cum cacteris sui ordinis coëfii- 
cientibus complere debeant unitatem ; sic D valerc diximus —1/30, 

• 1 1 2 7 2. . 1 ^ 

quia - H \ h - i+D) = 1. 

^ 9 2 3 15 9^ ^30 

Huius laterculi beneficio intra semi-quadrantem horae reperi, quôd potestates déc- 
ime sive quadrato-sursolida mille primorum numerorum ab unitate in summam 
collecta efiiciunt 

91 409 924 241 424 243 424 241 924 242 500." 



1.2. Définissons les nombres 6„ par la série génératrice 



1 - e-^ ^ n! 2 ^ 

On remarque que 

^ - 1 I V J^Q 

eS_i 2 (2p)! 

Voici les premiers valeurs : 
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1.3. On verra plus bas que Sr{n) est la valeur en a; = n du polynôme 



l<p<(r+l)/2 

Autrement dit, 



S'^{x) = B^{x) = ^ ( j^X"'', 5,(0) = 



§2. Une primitive et l'équation de Rota-BaLxter homogène 

2.1. Soit A une algèbre de fonctions j{x) "raisonnables", par exemple l'algèbre 
des polynômes M[a;] ou l'algèbre des fonctions dérivables. On désigne par D l'opé- 
rateur de dérivation sur et par 7 l'opérateur 

Af){x)= r mdt. 

Jo 

Il est clair que DI = id. En revanche, 

iD{f){x) = f{x)-m . 

2.2. Lemme. L'opérateur / satisfait l'équation 

/(/)J(ff) = I{I{f)g + fl{g)) . {RBH) 

Première preuve. Les dérivées des deux côtés coïncident puisque DI = id. De 
plus, les valeurs des deux côtés en sont 0, d'où l'assertion. 

Seconde preuve. Considérez l'intégrale de la fonction de deux variables f{t)g{u) 
sur le carre [0, x]"^ ; puis coupez ce carré en deux triangles. 

L'équation (RBH) sera appelée équation de Rota-Baxter homogène, cf. [Ro]. 
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§3. Une primitive discrète et l'équation de Rota - Baxter non-homogène 

3.1. Étant donnée une fonction / : N — > M définissons sa "primitive discrète" 
B{f) : N — > M par 

n 

B{f){n) = j2 m- 

i=l 

3.2. Lemme. L'opérateur B satisfait l'équation 

Bifg) + B{f)B{g) = B{B{f)g + fB{g)) . (RB) 

En effet, la valeur du premier membre de (RB) en n est 

n n n 

i=l i=l j=l 

Donc c'est un carré n x avec la diagonale doublée. D'autre part, le deuxième 
membre est 

n f i i N 

E E/œ^w+E/w^œ • 
i=i "-^=1 j=i j 

Il est aisé de voir que les deux expressions sont égales, en interprétant chaque terme 
Sj=i /(i)â'(^) + Sj=i /(^)â'(i) comme un chemin (de forme F) dans le même carré. 

3.3. Il est clair qu'on peut considérer l'anneau de polynômes R[x] comme un 
sous-anneau de l'anneau des applications / : N — > R. 

Lemme (Bernoulli). B{R[x]) C R[x]. 

En effet, on peut calculer les polynômes Sr{x) := B{x^) par récurrence, en 
utilisant (RB) : 

On a B{l){n) = n, donc B{1) — x. Ensuite, 

B{1 ■ 1) + B(1)B{1) = B(B{1)1 + IB{1)), 

i.e. 

x + x^ = 2B{x), 

d'où B{x) = (a;2 +a;)/2. 
Ensuite, 

B(l • x) + B{l)B{x) = B{B{l)x + lB{x)), 

i.e. 

+ = B{x^ + ) = 2^(3^ + ^) = 2^^^^ + ^^' 

d'où 
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cf. 1.1. Ainsi, les B(x*) pour i < r étant connus, on obtient B{x^~^^) en appliquant 
(RB) avec f = l,g = , ce qui prouve le lemme. 

3.4. Soit A une algèbre associative munie d'un opérateur linéaire B : A — > A 
vérifiant 

^iB{fg) + B{f)B{g) = B{B{f)g + fB{g)) (RB) 
(/Lt est un nombre). On introduit sur B une autre multiplication 

f*g = B{f)g + fB{g)-t,fg. 

En l'utilisant, on peut réécrire (RB) sous la forme équivalente 

B{f * g) = B{f)B{g) . {RB)' 

Autrement dit, B est un opérateur entrelaçant deux multiplications. 

3.5. Lemme. La multiplication * est associative. 
Se vérifie aisément à l'aide de (RB)'. 

3.6. Lemme. L'opérateur B satisfiait à (RB) pour * : 

^iB{f * g) + B{f) * B{g) = B{B{f) *g + f^ B{g)) . 

En effet, 

B{f * g) + B{f) * B{g) = B\f)B{g) + B{f)B^{g) . 

D'un autre côté, d'après {RB)' , on a B{B{f) * g)) = B^{f)B{g) et B{f * B{g)) = 
B{f)B^{g), d'où l'assertion. 
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§4. Formule sommatoire d'Euler - Maclaurin 



4.1. Soit A une algèbre commutative munie d'une dérivation D et d'un opérateur 
/ : A — > A ("une primitive") tel que DI = id^ et satisfaisant 

Hmg) = mf)9 + fiig)) • {RBH) 

On veut construire un opérateur B : A — > A ("une primitive discrète") de la 
forme 

B = 7(1 + /îi£) + /32£>^ + ...). /3ieM, (4.1.1) 

qui satisfait 

B{f9) + B{f)B{g) = B{B{f)g + fB{g)) . {RB) 
Essayons de trouver un par un les coefficients inconnus /3j, de l'équation (RB). 

4.2. On a à gauche : 

B{fg) = I{fg) + PiIDifg) + p2lD\fg) + ... 

= Hfg) + PiHDf ■g + f-Dg)+ p2l{D^f ■g + 2Df -Dg + f ■ D^g) + ... 
Ensuite, 

B{f)Big) = (/(/) + piIDif) + p2lD\f ) + ...)■ {lig) + 0^1 D{g) + p^ID^ (5) + • • • ) 

= I{f)I{g) + (3,iIDif)I{g) + I{f)ID{g)) 
+02{ID\f)I{g) + miD^g)) + (illD{f)ID{g) + ... 
= I{Iif)g + fIig)} 
+(3il{lD{f)g + Dfl{g) + I{f)Dg + ÎID(g)} 
+(32l{lD\f)g + D\f)I{g) + I{f)D^g + fID\g)} 
+fill{lD{f)D{g) + D{f)ID{g)} + . . . 

4.3. À droite : 

B{B{f)g + fB{g)) 

= l[l{î)9 + fl{9) + Pi{ID{f)g + f ■ ID{g)) + (i2{ID\j)g + / ■ ID^g)) + ... 
+/3i7£>|/(/)5 + fl{g) + 0i{ID{f)g + f ■ ID{g)) + . . . | 
Wi5'|/(/)5 + //(5) + •••} + •• • 
= l[l{f)9 + ÎI{9) + Pi{ID{f )9 + f ■ mg)) + p2{ID\f)g + f ■ ID\g)) + . . . 



+/3i7|/5 + I{f)D{g) + D{f)I{g) + fg^ 
+Ppl^D{f)g + ID{f)D{g) + D{f)ID{g) + fD{g)^ 
+(}2l{ D{f)g + 2fD{g) + I{f)D\g) + D\f)I{g) + 2D{f)g + fD{g) \ + 



(on a gardé les termes d'ordre < 2 en £)). 

4.4. En réalisant l'équation, les seuls ternies qui surviennent sont : I{fg), ce 
qui donne 

1 = 2/3i, (4.4.1) 

i.e. /?! = 1/2 = — 6i ; puis, 

I{D{f)g + fD{g)), qui donne 

A =3/32+/??, (4.4.2) 

d'où 

12 2-6 

De la même manière, les calculs à l'ordre 3 fournissent la valeur /îs = 0. Ils donnent 

à l'ordre 4, en regardant les termes I{D^(f)g + fD^{g)) (et aussi I{2D'^{f)D{g) + 
2D{f)D'^{g))), la rélation 

= 5/34 + /32, (4.4.3) 



d'où 



5 5-144 4! • 30 ■ 

Ceci nous amène à l'idée que 

I3„ = ^, n>2. 
n! 

Autrement dit, on attend que la série génératrice des nombres /3„ soit 

S 



1 - e-^ 

n=0 



4.5. Écrivons la série (4.1.1) comme 
où = ix{D) e IR[[£»]], ^(0) = 1. L'équation (RB) s'écrit alors : 

ifi{f)iK9) + iKfg) = iKiMf) -g + f- i^g)) ■ 

On pose : /x(/) = o, fi{g) =h,v = donc / = p{a),g = y{h). On obtient : 
Ia-Ih + In{iy{a)iy{b)) = In{Ia ■ v{b) + v{a)Ih) . 
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On multiplie les deiix côtés par i^D (en prenant en compte le fait que DI = id) : 
uialb + Ia-b) + i'{a)i'{b) = la ■ v{h) + v{a)Ib . 

Ensuite, on pose la — a, Ib = j3, donc a — Da, b ~ D(3 : 

vD{a(3) + yD{a)vD{P) = avD{l3) + uD{a)f3 . 

Il s'en suit que si l'on pose k = 1 — vD, alors 

K{aP) = K{a)K{f3) . (4.5.1) 

Il est clair que «(1) = 1. 

4.6. Maintenant supposons que notre algèbre A est l'anneau de polynômes 
A = R[x]. L'identité (4.5.1) implique que k est un automorphisme de A, donc il est 
de la forme k{x) = Ex + F. 

De plus, on impose une condition de normalisation 

B{l)=x, (4.6.1) 
d'où //^(l) = X, donc IJ.{1) — 1, donc J^(l) = 1, d'où 

k{x) = (1 - i^D){x) =x-l . 
Il vient que k = , d'où vD = 1 — e~^, u = {1 — e~^)D~^, donc 

D 



H{D) = 



1-e 



-D ' 



comme attendue. Donc 



4.7. Explicitement, on a 



B=^-^. (4.6.2) 



1-e-D 2 (2p)! 

p=i ^ 

Maintenant, pour f{x) € on a 

iD{f){x) = f{x)-m, 

donc 



m- 



On a montré que cet opérateur vérifie (RB) et (4.6.1), or, il n'existe qu'un seul 
opérateur de la sorte (l'unicité se voit tout de suite par récurrence), celui qui à 
/ e R[x] fait correspondre B{f) e R[x] tel que B{f){n) = ElLi /(«)• 
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Il vient le 

4.8. Théorème (Euler - Maclaurin). Pour tout / G R[x] et n e N on a 

/(l) + /(2) + ... + /(n) 

En l'appliquant à f{x) = , on obtient 

V + 2"- + ... + rf = Sr{n), 
où le polynôme Sr{x) est défini par 

l<p<(r+l)/2 ^ ^ ^ ^ 

4.9. Appliquons (RB) à / = 1, g = : 

B{x'') + xB{x'') = B{x''+'^ + B{x'')), 

i.e. 

B{x''+'^) = {x + l)B{x'') - B{B{x'')) . 
Cette identité est équivalente à une identité classique pour les nombres de Bernoulli : 

/'2n\ 

(2n + l)b2n = ( „ ) b2pb2n-2p, 



p=l ^ ' 

cf. [Bo], Ch. VI, §2, Exercice 2) ; [R] (c'est le premier article publié de Ramanujan). 
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Deuxième Partie 
UNE FORMULE DE RAMANUJAN 
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§1. Fonction rj de Dedekind 



1.1. Il semble que Riemann ait lu assez attentivement les Fundamenta de Ja- 
cobi. Dans les papiers de Riemann, on a trouvé un Additamentum ad §"™ 40 de 
"Fundamenta", [R]. Richard Dedekind a écrit un commentaire sur ces fragments, 
[D], où il introduit la fonction //(r) : 

oo 

,7(t) = -Q (i_^n)^ g = e2--, (1.1.1) 

n=l 

oii \q\ < 1, i.e. 9t > 0, et étudie sa loi de transformation par rapport aux trans- 
formations de Moebius r i— > (ar + b)/{cT + d). Le théorème suivant en est un cas 
particulier. 

1.2. Théorème. La fonction ?7(t) satisfait l'équation 

r?(-l/T) = VWivir) . (1.2.1) 

Démonstration^ d'après Cari Ludwig Siegel, [S]. En prenant le logarithme naturel, 

oo oo nm i 

![!i-iogr?(T)=-y iog(i-g")= y —=y , • 

^2 ^ ' ^ 6V y ; m m(q-"^-l) 

n=l n,m=l m=l ^ ' 

Prenons le logarithme de (1-2.1) : 

log7?(-l/r) = ilog(T/z) +l0g77(T), 

ou 

■ïïiT , , , TTZT"-^ , / -, , N 1, / ,.x 7ri(r + r~-^) 
— - logr;(r) = log77(-l/r) + - log(r/i) + . 

Donc (1.2.1) est équivalente à : 

1, , 7ri(r + T-i) 1 ( 1 1 \ 

m=l ^ ' 
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1.3. Une fonction intéressante : cotz. On pose y = e". Alors : 

_ cos^ _ {y-'^ +y)/^ _ ■ y^^ + y _ . y + y^^ 
sin^ {y-^-y)/2i y~^-y y-y~^ 



y-2-l)-*-(^+y2_i 
Donc limy_>o cot z = —i et linij^^oo cot z = i. De là : 



lim cot((n + l/2)z) = si 3z > (1.3.2a) 

n— J-oo 

et 

lim cot((n + l/2)z) = i si 9^ < . (1.3.26) 

n—^oo 

1.4. On pose f{z) = cot z cot z/t et on considère la fonction gn{z) = z~^f{vz) 
ovlU = [n + l/2)7r, n = 0, 1, . . . Soit C le contour du parallélogramme de sommets 

l,r,-l,-r. 

Quels sont les pôles de ^«(z)? On a : 

cosi^z cosvz/t 

9n{z) = —. : 1- . 

zsiauz smvz/T 

Donc on a : 

(a) des pôles simples en ^; = ^irm/v, m = 1,2,..., avec les résidus 

cot(7rm/r) 



reSjg— -j-TT-m/i^ 9n{z) 



nm 



(b) des pôles simples en ^; = ±.-KmT/v, m — 1,2,..., avec les résidus 

, . cot(7rmT) 

reSjj^-i-TrTOT/i' 9n\z) — • 

(c) Enfin, en 2; = on a : 

1 J_ T_ 1 - iy^z^/2 + ... 1 - iy^zy2T^ + ... 
^"■^^^ ~ z' vz' vz' \- v-^z^l^ + . . . ■ 1 - z/22V6r2 + . . . 

v^z"^ \ v^z^ 



" •(l-^-(l + r-) + ... 



d'où 
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Par la formule des résidus de Cauchy, 

dz T + T-'^ 2^1 



1 f ^, ,dz r + T-i 2 ^ 1 , 
/ .f('^z) — = 1 — > — (cotTrmr + cotTrm/r) . 

^ m=l 

On remarque que 

cot Trmr + cot nm/T = —2i j 
cf. (1.3.1), d'où 



1 1 

-27Timr ^ g27rim/r 1 / ' 



"-^ m=l ^ ' 



1.5. Maintenant faisons tendre n à l'infini dans (1.4.1). Soit £i = {tsz = 0} et 
£2 la droite qui passe par et t. D'après (1.3.2a,b), 

lim„^oo cot = —i si z est au-dessus de £1 ; lim„^oo cot vz = i si z est au- 
dessous de ^1 et 

lim„^oo cot z^2;/r = i si ^; est à droite de £2 ; lim„^oo cot î/2;/r = — i si 2; est à 
gauche de £2- 

11 s'en suit que sur le côté (1,t) de C (sans les sommets) la valeur limite 

lim„^2c cot vz cot i>z/t = —i-i= 1. 

De même, sur les côtés (r, — 1), (— 1,— r) et (— r, 1) les valeurs limites sont 
-1,1,-1. 

De là, 

"^°<=Jc z \Ji J_i J_^J z 

= logr — 7r + logT + log(— r) — TT — 27r + log(— r) = 41ogT — 27r = 41og(T/z) . (1.5.1) 
Donc en passant à la limite n ^ 00 dans (1.4.1), on obtient : 

41og(TA) + =S}_^ - l _ ^ - ^^^.^^^ _ 1 . 

m=l ^ ' 

En divisant par 8, on obtient la formule cherchée (1.2.2), QED. 
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§2. Une formule de Schlômilch 



2.1. Théorème, [Sch], [Ram]. 

y , " =1-^. (2.1.1) 

2.2. Démonstration de Srinivasa Ramanujan, [Ram], (18), p. 32. On prend 
T = ia dans (1.2.1), où a est un nombre réel, a > : 

•OC' OO 

g-,r/12a (1 _ 6-2^^"/") = . g-^/^^ Il ~ ^-^^na^ ^ 

n=l n=l 

En prenant le logarithme, 

OO 



12a ^ °' ' 2 12 

n=l n=l 



y: log(l - e-W») = 1^ _ ^ + ^ iog(i - e-2-«) . 

En prenant la dérivée, 

TT ^ (27rn/a2) ■ g-^W" _ J_ _ ^ ^ 27rne-2'^"/» 
1202 ~ ^ 1 - e-2'r"/« ~ 2â ~ 12 ^ T 

n=l n— 1 



OU bien 

.IL(Q-2 + i)_l=27rV n- f +^ V (2.2.1) 



12 '2a ^ \(p.-^nla_y g27™o _ X ; 

n=l ^ ' 

Sous une forme plus symétrique, 

!L(^l±^_i = 2.y f ^ + , . (2.2.2) 

]^2 2 / ^ y g27m/a _ ]^ g27ma — 1 / 



En posant a = 1, on arrive à (2.1.1). 
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§3. Développements eulériens de sin et de cot 

3.1. On suit Bourbaki, [B], Chapitre VI, §2. 
Lemme. On a pour n e Z, n > : 

n-l 

sinn^ = 2"""^ sin(^ + kn/n) . 
fe=o 



En effet, 



1 -mz 

sinnz = -(e'"^ - e"^"^) ^ (e^™^ - 1) 



2i J--»-' ' 2i „ 

piz+Tvip/n p—iz—'ïïipin 

= (2z)"-i n e-^'V" ]J f _i 

p=0 p=0 



Or, 



d'oii l'assertion. 

3.2. En divisant par sinz et en faisant tendre z vers 0, on obtient 

n-l 

JJ sin(î37r/n) = «2^-" . 

3.3. Soit n = 2m + 1 impair. On a : sin(n(2: + 7r/2)) = sin(n2; + 7r/2 + m7r) = 
(— 1)™ cosnz, d'oîi, en remplaçant z par z + 7r/2 dans 3.1, 

n-l 

cosn^ = (-l)'"2"-i JJ cos(^ +p7r/n), 

p=0 

donc 

n-l 

cotn^; = (-l)™2"-i JJ cot(^; +p7r/n), 

p=0 

que l'on peut réécrire comme 

m 

cotn^= (-l)™2"-i cot(0 - p7r/n) . (3.3.1) 
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3.4. On a : 

cos(a + 6) cos a cos 6 — sin a sin 6 



cot(a + b) 



sin(a + 6) sin a cos 6 + cos a sin 6 

1 — tan a tan h 



tan a + tan b 

Donc 

tan2; — tan(p7r/n) 

p=—m " ' ' 

Ceci est une fraction rationelle dont le numérateur est de degré n — 1 en m = tan z 
et le dénominateur est de degré n, ayant les racines simples. Il s'en suit qu'on peut 
écrire une décomposition en éléments simples : 



cot nz 



tn 

11. — tar\( 



u — tan(»7r/n) 

p=—m ^ ' ' 



avec 



ttp = lim cotn^ • (tan^ — tan(p7r/n)) 

>p7r/n 

cosn^; sin(2; — n7r/n) 

— lim 



z^p-iv/n sinn^; cos2;cos(p7r/n) 

1 ,. cos(nft, + p7r) sin/i 1 {—ly sinh 

— Ô7 TT 7 ^ — = — Ô7 TT ^1™^ 7^ — r 

cos^(p7r/n) fc^o sm{nh + p-ïï) cos-' [pw / n) h^o [—1 )p sin nh 

1 



ncos^(p7r/n) 

Donc 

m ^ 

cot nz = 



n cos^ (jm/ n) (tan z — tan(p7r/ n)) 



En remplaçant z par z/n, 

m 

cot z = 



ncos^(f)7r/n)(tan(^;/n) — tan(p7r/n)) 



1 ^ 1_ 



2ta.n{z/n) 



ntan{z/n) n cos^ {jm /n) t&v? [z / n) — tax)? [p-K / n) 

1 2ntan(z/n) 



E 



ntan(z/n) cos2(p7r/n)(ntan(z/n))^ — (nsin(j37r/n))2 



p. 

On a donc démontré le 



3.5. Théorème. Pour tout n = 2m + 1 impair 

1 2ntan(z/n) 



cot z = — — — -— + 

n. tan l z n.) ' 



ntan(^/n) ^ cos2(pVn)(ntan(^/n))2 - (nsin(p7r/n))2 " 
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En faisant m — > oo, on arrive à : 
3.6. Théorème. 

2z 



1 ^ 

cot z = — h 



z ^-^ z"^ — »^7r^ 
p=i 

3.7. Revenons au développement de sinus 3.1. Supposons toujours que n 
2m + 1 est impair. Alors 3.1 peut s'écrire 

m 

sinn^ = (-l)"2"-i sin(z-p7r/n) 

p=— m 

m 

= (-1)™2"'~^ sin z sin(2; - pn/n) sm.{z + piï/n) . 

On vérifie aisément la formule suivante : 

sin^(a + b) — sin^(a — b) = sin 2a sin 26, 

d'où 

sinosin6 = sin^((a + h)/2) - sin^((a - 6)/2) . 

Il s'en suit, 

sin(z — pir/n) sm.{z + pn/n) = sin^ z — sin^ {pn /n) , 

d'où 

m 

sin 712; = 2"~^ sinz (sin^(p7r/n) — sin^ z) . 
p=i 

Or, d'après 3.2, 

m 

d'où 

m 

sinn^ = nsin2; (1 ^ (sin^ z/ sin^(p7r/n))) . 

p=i 

En remplaçant z par z/n, on arrive au 

3.8 Théorème. Si n = 2m + 1 est impair alors 



sm^; = nsin(2:/n) I I 11 ^ 

V sin [kw 



/n) 



Maintenant si l'on fait tendre m vers l'infini, on obtient le 
3.9. Théorème. 

oo ^ 

sin = 2; • I 1 



p=i 



p'^ir^ 
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(Convergence uniforme dans des sous-ensembles compacts.) 
Application aux nombres de Bernoulli 
3.10. On a : 

. e-^ + . e^ + l 
cotiiz 2) = i ■ = —i ■ , 

^ ' ^ e-^-e^ e^-l 

d'où : 

z z ( _ + 1 \ z iz , . 

e^-1 2V e^-iy 2 2 ^ ' ^ 
On rappelle que les nombres de Bernoulli sont définis par : 



(2n)! 

71=1 ^ ' 



3.11. Le développement de cot nous dit : 

1 ^ 2z 

Maintenant : 



n=l 



20 2^ 1 2z 



' fe=0 



E 



^2fe 



-2E 



k=\ 



\z\ < n). En échangeant l'ordre de sommations, il s'en suit : 

4. l oV^ "^2*; 2fe-l 

cotz=--2}_^ -^z , 



ou 



Donc 



oo ^ 

= E ~k ■ 

n=l 



k=l 



1 _ £ + V f-n^-i z 



/c=l 

3.12. En comparant avec 3.10, 
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OU 

\2n 



n> 1. 



52„ = (-ir-^|^62„, 



§4. Une formule de Ramanujan 



4.1. On agit à la Eisenstein. On suit [A], Chapitre II, no. 10. Commençons par 
le développement de cot : 



TT cot nu 



11. ^-^ 



u ^ — ' , \ M + m m 



1 1 



= i V { ^ — 

7/ \ 11. -\- m 01 — 

On pose w = e^^*^ ; alors 



u \u-\-m u — m 

m=l ^ 



cot nu = i = i ■ 1 H \ = —i + 2i y w , 

w — 1 \ w — 1 J 



71=1 



si \w\ < 1, i.e. ?su > 0. Il s'en suit, 

1 ^ 1 1 ^ ■ n n 

- + > 1 = -m - 2m > li;" . 

u ^-—'^ \u + m u — m I 

m— 1 ^ 



4.2. On dérive p fois par rapport à u ; puisque {d/duY{w) = {2ni)Pw, on a : 

oo oo 

{-ifp\ y - — =-{2ni)p+^y 

ra=-oo ^ ' k=l 



On pose u = nr, n > 0, Sr > 0, 



oo 

2kn'irir 



m= — oo ^ ' K=\ 

et l'on réalise la somme sur n : 

°° oo 1 oo 2kniT 

n=l m=— oo ^ ^ fe=l 

(attention : on a changé l'ordre des sommations à droite.) 
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4.3. Maintenant supposons que p = 21 — 1 est impair et p > 2 (i.e. l > 2). On 
peut alors réécrire (4.2.1) : 



2 ^ (m + nryi ~ ^ ~ (2Z - 1)! ^ T 
On utilise la notation 



Im . 4- r 



(m + nr)*^ 

pour les séries d'Eisenstein. 

4.4. Considérons le cas spécial t = i : 



1)! ^ 



(2/ - 1)! e^fe'^ - 1 ■ 

On rappelle en revanche que 

cf. 3.12. Il s'en suit : 

^2^^ - (-1) + ^ • 

4.5. Supposons que l = 2j + 1 est impair. Alors 



(m + ni)^^ (— mi + n) 

donc E2i{i) = 0. Il découle que 

fc=i 

dans ce cas. Ceci est une formule de Ramanujan. 

4.6. En général, on définit les fonctions de Weierstrass : 

a(u)=a(a;i,a;2;u)=un ' (|l - ^^e"/'^+"' 

oii w = mu)\ + mo^ et 

n'= n 

(m,n)eZ2-{(0,0)} 



1 u 

- + ^ 



Cette fonction est analogue de sinu. Ensuite, 

c(-)=4t=E'(— 

^ ^ a{u) ^ \u- 
analogue de cot u ; et 

n.) = -<'(») = E'{ïï^-ii}. 

analogue de — cosec^w. On a alors le développement de Laurent en : 

au) = -- E4U^ - Equ^ - Esu^ - ... , 
u 

oh ^ 

En = En{uJi,U-2) = y] ' — • 

Donc 

Viu) = \+ ^Eiu'^ + ^Equ^ + 7Es,u^ + ... 
La fonction V{u) satisfait les équations différentielles 

V'^{u)=AV\u)-g2V{u)-gs , 
V"{u) = 6V^{u)-gs/2, 

où. 

52 = 60Ei, 53 = UOEq . 

4.7. Le cas du réseau Gaussien {oji,co2) = a été traité par Hurwitz, 
On considère la fonction de Weierstrass qui satisfait l'équation différentielle 

= 4P^(m) - 4:V{u), 

donc (/2 = 1, 53 = 0. On introduit la période correspondante : 

„ / dx 
u} = 2 

cf. une définition de tt : 



7r = 2 





1 dx 



On définit alors les nombres rationcls £"„ par 

1 2'^Ei 2^E2 2'*"£;„ 



riu) = 



8 6! 4n (4n-2)! 

On a, El = 1/10 et satisfait une relation de récurrence 

(2n-3)(16n^-l) g (^^ - - 4fc - 1) ■ 
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Alors 



E' 



4n 



1 _ (2w) 
(r + iâ-)*" ~ (4n)! 



§5. Une intégrale de Legendre 



5.1. On rappelle que 



t t t t t b2ni'^^ 



On peut donc poser 61 = —1/2. 
5.2. Théorème (Legendre). 



°° sinax _ 1 e° + 1 _ J_ 
e^»^»; _i^~4e«-l 2a' 



On donne deux démonstrations. 

5.3. La première démonstration utilise le développement de cotangent. 
On a : 



(a; > 0), d'oii 



g27rx ^ 

n=0 n=l 



/.oo sin /'OO 
Jo e2-- - 1 



Or, 



ou 



/•oo -1 /"OC 

/ sinaa: e-2'^"^da;= - / (( 



^iax—27Tnx ^ ^iax—27Tnx 



Iq ia — 2-Kn 

Donc 



1 27rn + ia 



2-ïïn — ia + An'^n'^ 



i 



00 

sinaa; e~^''"''c?a; = 



d'où 



00 

a 
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Rappelons : 



De là 



Or, 



donc 



2a 1 
/ ~^ = cot a . 

a-^ — 'K'^n'' a 

n=l 



oo 



n=l n=l ' ra=l ^ ' ' 

= -l(cot(^a/2)-A)=-lcot(za/2)-^. 



_ cos(m/2) _ i(e-°/^ + e°/^) _ z(l + e°) 
cot(^a/2) - ^r^^^-y^ - g_„/2 _ ' i _ g» ' 



■ cot(za/2) = 



4z ^ ' ' 4 e« - 1 
quod erat demonstrandum. 



5.4. La deuxième démonstration utilise la formule de Cauchy ; elle a été proposée 
comme un exercice dans [WW], Ch. 6, 6.4, Example 2. Le calcul a été fait par 
Nabil Rachdi. 

On définit : 
Alors 

c^o 2i 

Considérons le contour "rectangulaire" F = r(e,iî) suivant : 
6 

r = [J Ti = {e< z < R}U {z = R + it\ <t <1}U {z = t + i\ R>t> e}U 

i=l 

U{z = i + ee'> >e> -tt/2} U {z = it\ 1 - e > f > e} U {z = ee*^| 7r/2 > é» > 0} . 
On pose : 

= ■ 

Puisque e^^^^ = 1 ssi z = ni, n e Z, cette fonction n'a pas de singularités à 
l'intérieur de F, donc 



= / f{z)dz = V / f{z)dz 
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Calculons les intégrales Jp f{z)dz séparément. On a : 

/ f{z)dz = Iie,R) . 

De même, 



f{z)dz = - ^-^^^-^dx = -e-'^I.{e,R). 



Ensuite, 

f{z)dz= I 

ir2{R) 

Les intégrales sur les quarts de cercles : 

-V2 ia(î+ee''') 



/ f{z)dz = / quand iî ^ oo . 



/ f{z)dz= / , iee^^d^ 



Or, la fonction sous l'intégrale 



e27ree*e _ 1 2-K 

d'où 



quand e ^ 0, 



/r4(e) 

De même, pour Jp on trouve : 



/ f{z)dz — > " quand e ^ . 

Jr.U) 4 



/ f{z)dz — > — quand e — > . 

■/r6(e) 4 



5.5. Finalement, il reste à traiter l'intégrale Jp^. On a : 

fl—i ^ia-it pl — e ^—at 

f{z)dz = - I 
Par la formule de Cauchy, 



= (1 - e-)/.!») - j(l + e-) + J„(«) + 



d'où, en posant y = e". 



25 



Il s'en suit : 



Or, 



{x = -t+ 1) 



-a(e) = -i^ 



1-e gat 



1-6 6 



• {—dx) = —iy 



l — e ^'Kîx ^ g— aa; 

1-e 



dx 



-ly 



e "^rfa; . 



La dernière intégrale 

f-l-e 



e-«^rfa;= / e-"" dx + o(e) = - - 

lo ^ ' a 



^ / N 1 — y ^ , , 
+ o(e) = f — + o(e) 





a 



Ainsi, 



J_„(e) = -yJ,(e) + '-^^-^ + o{e) . 



(5.5.1) 



On obtient : 



7„(e)-/_„(e) = - 



l + y'^ _ 1 + 2/ 



1-y 



Ici : 



et 



_ Ja(^) = _ 1 _ ^ + = _1 + . (5.5.2) 

1 — y 1 — 1 — y y~i û 

On peut voir en (5.5.1), ou en la formule équivalente (5.5.2), une équation fonc- 
tionelle pour la fonction Ja(e) ; remarquons que l'intégrale Ja(e) diverge quand 
e ^ 0. En revenant k la, on obtient : 

4(e)-/-«(e) = --L__+,(e), 



d'oii 



Ue)-I-aie) 1 y + l 1 . . 



2z 4 y - 1 2a 

En faisant tendre e vers zéro, on obtient la valeur de l'intégrale de Legendre. 



5.6. Théorème. Pour n > 1, 



oo j.2n— 1 

■n-li,_ _ A„ I Jj. 

p27rt _ 1 • 



(5.6.1) 
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On peut considérer cela comme une deuxième définition des nombres de BernouUi 
(Jacob BernouUi, Ars conjecfMudi, 1713, p. 97). 

On en donne deux démonstrations. 

La première démonstration utilise les valeurs de C,{s) en points positifs pairs 
(donc le développement de cot) : 




{x = 2-Kkt) 

= V / (a;/(27rfc))2"-ie-^rfa;/(27rA;) 
= (27r)-"r(2n) Y. ^ = (2'r)-"(2n - ly.S^n 

(voir 3.12) 

^ ' An ' 



5.7. La deuxième démonstration utilise l'intégrale de Legendre (avec la preuve 
par la formule de Cauchy), cf. [WW], 7.2 : 



°° sinax , 1 i 1 ^ , i^aY 

, „ V^j,^ 



-dx = —- — h - cotza : 



e'^^-l"" 2a ' 2 2a ^ ^""(2n)! 

" 71=1 ^ ' 

En dérivant 2n fois et en posant a = et x = 2t, on en déduit (5.6.1). 
En particulier, si n est impair, 

b2n r f^-^ 



Jo 



An Jo e^'^* - 1 

cf. 4.5. On arrivé ainsi à l'assertion : 

5.8. Theorema pulcherissimum. Si n > 1 est un entier impair, alors 



_ _2n 

g27rfe _ 1 ~ 
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